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ABSTRACT 

We study the dependence of the surface magnetic fields of radio pulsars on the choice 
of Equations of State, pulsar masses and the values of the angle between the magnetic axis 
and the spin axis of the pulsars within simple dipole model. We show that the values of the 
surface magnetic field can be even order of magnitude different from its canonical values. 
This difference will effect any magnetosphere related model to explain observational features 
of radio pulsars and magnetars. We find a significant difference of the value of the surface 
magnetic field from the commonly quoted value for the faster member of the double pulsar 
system, i.e. PSR J0737-3039A as here both the mass of the pulsar and the angle between the 
magnetic axis and the spin axis are known. Our study reveals the importance of constraining 
the dense matter Equations of State in pulsar astrophysics as well as hints an alternative way 
to constrain these by independent determination of the pulsar magnetic field. 
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1 INTRODUCTION 

Magnetized rotating neutron stars can be categorized into different 
categories based on their emission properties, like AXPs, SGRs, 
radio pulsars, RRATS etc. Although almost every day some new 
interesting discoveries relating these objects are being made using 
the advanced technologies, it is unfortunate that the true description 
of the dense matter constituting these objects is still not well under- 
stood. A number of Equations of State (EsoS) are available in the 
literature trying to describe the state of the matter in such extreme 
condition as inside neutron stars, see Lattimer & Prakash (2007) 
for a few examples of the EsoS. Better knowledge about the dense 
matter EsoS will in the long run help us to understand observational 
features of neutron stars in a better extent. 

Presently there are a number of approaches trying to constrain 
the dense matter EsoS through astronomical observations of com- 
pact stars. The usual approach is to determine the mass and the 
radius of the stars with the help of various observational features 
like gravitational redshifts (z) from spectral lines, cooling charac- 
teristics, kHz quasi-periodic oscillations (QPO) etc. (Lattimer & 
Prakash 2007, Li et al. 1999, Ozel 2006, Zhang et al. 2007). But 
these methods are not foolproof, e.g. the value of z used in Ozel's 
analysis of EXO 0748—676 can not be reproduced as mentioned by 
Klahn et al. (2006). Moreover, to constrain EsoS from QPO obser- 
vations, one need to believe in some specific model of QPO which 
is again a subject of debate. Another alternative method might be 



the measurement of the moment of inertia from the faster compo- 
nent (A) of the double pulsar system PSR J0737-3039A/B (Lat- 
timer & Schutz 2005, Bagchi et al. 2009). As some high mass 
stars like PSR J1903+0327, EXO 0748-676 etc. prefer neutron 
star models and some other stars like 4U 1728-34 (Li et al. 1999), 
prefer strange star models, it is possible that both family coexist 
in nature. But even then, we need some constrains as there are a 
number of EsoS for neutron stars and also for strange stars. 

Because of the lack of any strong constrains on the dense mat- 
ter EsoS till date, it is interesting to study the effect different choice 
of EsoS on the stellar properties. This is our motivation of studying 
the effect of EsoS in the determination of the values of the pulsar 
magnetic field and consequent results. Here we have chosen five 
EsoS, three for the neutron matter and two for the strange quark 
matter. We discuss the model of the pulsar magnetic field in section 
|2]in short and then display our results in section[3] Discussions and 
conclusions are given in section|4]and section|5]respectively. 



2 MAGNETIC FIELD 

For a radio pulsar, the value of the surface magnetic field (Bs) is 
estimated by equating the spin down luminosity with the dipole 
radiation power and one gets the following expression : 
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where / is the moment of inertia of the pulsar, R is its radius, P is 
the spin period, P is the rate of change of the spin period, a is the 
angle between the magnetic axis and the spin axis of the pulsar and 

c is the speed of light. Bjaa is defined as ^J^^^^- 

Due to the lack of knowledge about the actual value of 
a, it is usually assumed to be 90 °. Moreover, / is taken 
as 10 gm cm , i? as 10 km. Then -B/ac becomes 3.2 x 
'^'^^sec"^''^. So Eqn. ([TJ can be rewritten as : 
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Eqn. ^ is commonly used to calculate by putting the values of 
P in seconds and P as dimensionless. 

The problem with the above simplification is that, for a fixed 
value of the pulsar mass (M), one gets widely different values of 
R and I by using different EsoS for the matter. Moreover all the 
pulsars do not have the same value of M. In short, 7? and I can 
be different from their canonical values 10 km and 10*^ gm cm^ 
respectively depending upon M and the choice of the EoS. a can 
also have any value between ° and 90 °. As a result the value 
of Bfac will be different from 3.2 x lO^'' gm^/^cm"^/^sec"^/2. 
For this reason, it is worthwhile to study the dependence of the 
value of Bs on M, a and EsoS. In the next section, we study this 
dependence using Eqn. ([TJ. 



3 RESULTS 

Among numerous EsoS for the dense matter, we have chosen five - 
three for the neutron matter (namely EoS BPALI2, EoS APR and 
EoS MET 17 which are also used by Bejger et al. 2005) and two 
for the strange quark matter. One in the latter group is is EoS A 
(SSA) of Bagchi et al. (2006) and the other one is from Bag model 
(BAGl) with model parameters as : B = 60.0 MeV/fm^, lUs = 
150.0 MeV, rriu = md = 0, = 0.17 where B is the Bag pa- 
rameter, rua, rUu, TUd are masses of s, u and d quarks respectively. 
These five EsoS are of widely different stiffness and so they give 
significantly different plots in the M — R plane. In Fig.[T] we show 
the M — R and Al — I plots for all of these EsoS. It is worthwhile 
to mention here that MET 1 7 strongly overestimate and BPAL12 
strongly underestimate the nuclear matter incompressibility. APR 
seems to be a better description of the nuclear matter. 

We have calculated the moment of inertia using the formalism 
derived by Kalogera et al. (1999) within Hartle's approximations. 
We have found that even for an angular rotational frequency (fi) as 
high as 5000 sec~^, Bfac obtained by this method differs maxi- 
mum upto ~ 10% from that obtained by using proper codes for ro- 
tating stars (RNS code, written by Nikolaos Stergioulas and avail- 
able at http://www.gravity.phys.uwm.edu/rns/) for all of the EsoS 
mentioned earlier; the lower is the value of Q,, the smaller is the 
difference. So for radio pulsars, the use of Hartle's approximations 
is well justified as the fastest pulsar yet discovered is PSR J1748 
- 2446ad with ^ = 4500 sec"^ (see ATNF pulsar catalog at 
|http ://www.atnf. csiro.au/research/pulsar/psrcat/). 

We have plotted the variation of Bfac with pulsar masses for 
different EsoS (Eigs.[2]and[3l( and different values of a (1°, 30°, 
60° and 90°). We find that for the strange quark matter EsoS, the 
value of Bfac is usually large. Eor EoS SSA, the value of Bfac is 
greater than 3.2 X 10" for any values of a and M. For EoS BAG I, 
Bfac is greater than 3.2 x lO" for a < 90°, at a = 90°, Bfac ^ 
3.2 x 10^®. Eor the neutron matter EsoS, -B/ac is greater than 3.2 x 
10" for very low values of a, e.g. a = 1° (for any values of M). 
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Figure 1. Panel (a) shows the mass - radius curves for vaiious EsoS used in 
this work. Solid lines ai'e for the strange quai'k matter EsoS whereas dashed 
lines are for the neutron matter EsoS. Panel (b) shows the moment of inertia 
- radius curves for the same EsoS. 



But for a > 30°, -B/ac can be smaller than 3.2 x 10 depending 
upon the values of M. For EoS BPAL12, at a = 30°, Bfac < 
3.2 X 10" if M < O.9M0, at a = 60°, Bfac < 3.2 x lO" if 
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< 3.2 X 10" if A'/ < 1.4M( 



M < 1.3Mq, at q = 90° 
Eor EoS APR, at a = 30°, Bfac > 3.2 x 10" for any values of 
M, ata = 60°, Bfac < 3.2 x lO" if M < 1.6Mq, at a = 90°, 
Bfac < 3.2 X 10" if Af < 1.8Af0.EorEoSMET17, ata = 30°, 
Bfac < 3.2x 10" if A// < 1.5A/Q,butfora = 60° anda = 90°, 
Bfac < 3.2 X 10" for any values of M. 

For a fixed M and fixed EoS, Bfac is nunimum (Bfac. min) 
for a — 90°, giving the minimum possible value of the Bs 
(Bs, min) for that Af and that EoS. Erom Eigs.|2]and[3] it is clear 



that BfS^, 
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Let A-Bs|a/ denote the change in -Bs for a change Aa in a 
keeping M fixed and ABs\ct denote the change in Bs for a change 
A-Af in M keeping a fixed. Then for M being in the range of 

0.8 - 1.5 Mq, ABsir-^^" > ABslfs^ > AB,|^^« > 
ABsl'J^^'' > AS, I 
EoS; ABsjff^ > AS. 
AB. " 

Our study clearly reveals that in addition to constraining the 
EoS, the knowledge of a is also necessary to know the actual 
value of Bs for a pulsar of known A/, P and P. Sometimes it is 
possible to determine a observationally. As an example, for PSR 
J0737-3039A, Jenet & Ransom (2004) have found a to be either 
1.6°±1.3°orl4°±2°by careful study of the pulse profile while 
polarimetric study by Demorest et al. (2004) gives a = 4 ° ±3 °. 
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Figure 2. Variation of B/ac with M and a for different strange quark mat- 
ter EsoS. The unit of B/ac is gm^/-^cm~^''^sec~'^/^. 
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Table 1. Surface magnetic field (Bs) for PSR J0737-3039A with different 
values of a as estimated by Jenet & Ransom (2004) and Demorest et al. 
(2004). The canonical value of Bs for this pulsar is 0.064 X lO" G (us- 
ing Bfac = 3.2 X 10^^ gm'^/^cm^^/^sec"^/^). Pulsar parameters are 
M = 1.3381 Mq, P = 22.699 ms and P = 1.7599 x W'^^, 

xTpP = 1.9987 X 10-1" scci/2. 



Ct= 1.0 degrees. 
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EoS J, \ a ^ 


1.6° 


4.0° 


14.0° 


SSA 


4.577 


1.829 


0.528 


BAGl 


2.378 


0.955 


0.275 


BPAL12 


1.999 


0.803 


0.232 


APR 


1.763 


0.706 


0.204 


MFT17 


1.099 


0.440 


0.126 



In tabled we present the values of Bs for PSR J0737-3039A ob- 
tained from Eqn. ^ with different EsoS. We find that Bs can be 
even orders of magnitude higher than the canonical value obtained 
by using Eqn. 



4 DISCUSSIONS 

We have seen that the actual value of Bs can vary from its canonical 
value even orders of magnitude depending on the EoS, AI and a. 
This difference may effect various calculations involving the value 
of Bs like the study of Alfven QPOs in magnetars (Sotani, Kokko- 
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Figure 3. Variation of S/ac with M and a for different neutron matter 
EsoS. The unit of Bj-^^ is gm^/^cm— ^/-^sec— ^Z'^. 



tas & Stergioulas 2008), the study of timing properties of magne- 
tars (Ozel 2002), the study of magnetorotational collapse (Lipunov 
Gorbovskoy 2008) and various observed parameters of magnetars 
(Malov & Machabeli 2006) and so on. 

In Table [T] we showed the possible departure of Bs of 
the faster member (A) of the double pulsar system PSR J0737- 
3039A/B from the canonical value. Presently we can not do the 
same for the slower member PSR B as for this source, the value of 
a is unknown. But from Fig[2]and[3] we can anticipate the value of 
Bs to be significantly different from its canonical value. This fact 
will effect our understanding of the observed features of the dou- 
ble pulsar system in a great deal. As an example, once in a orbit, 
PSR A gets eclipsed by the magnetosphere of PSR B (which has a 
larger magnetosphere due to the slower rate of rotation). Lyutikov 
& Thompson (2005) modeled the eclipse as the attenuation of PSR 
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A's radio beams by the synchrotron absorption in the closed mag- 
netic field lines of PSR B. This model gives absorption coefficients 
as a function of the magnetic field of the pulsar B. Intensity and 
the polarization properties of the radio beam of PSR A during the 
eclipse is a function of the absorption coefficients and thus a func- 
tion of the magnetic field of PSR B. So any study of the eclipse 
light curve or the polarization properties to extract pulsar parame- 
ters might be effected by the use of the canonical value of Bs of 
PSRB. 

There is a common practice of drawing (\?o-magnetic field 
lines in the PP diagram for the pulsars using B/ac = 3.2 x 10^^. 
These line should be considered as iso-Bs, min lines (as a — 90°), 
instead of iso-Bs lines. Even then, only EoS BAGl gives these 
values of Bs, min for M = 0.8 — 1.4Mq as here B/ac, mm — 
3.2 X 10^^ EoS SSA will give higher values; EsoS BPAL12, APR 
and MFT17 will give lower values in this mass range. Following the 
same logic, the lowest value of Bs. min observed so far (as men- 
tioned in Ray Mandal, 2006; see their Fig 1) will shift to a higher 
value for EoS SSA and to some lower values for EsoS BPAL12, 
APR and MFT17 with M = 0.8 - IAMq. 

Eqn. ^ is derived by equating the pulsar spin down radiation 
with the magnetic dipole radiation power. This model predicts the 
pulsar breaking index n = 0,0, /Q? — 3. For some pulsars, the 
values of n have been measured observationally and most of these 
are found to be less than 3 e.g. Crab pulsar has n — 2.515 ± 0.005, 
PSR B 1509-58 has n = 2.837 ± 0.001 and PSR B0540-69 has 
n — 1.81 ± 0.07 (Lyne & Smith 2006 and references there in). 
This fact hints to the necessity of an improved model for the spin 
down of pulsars. Xu & Qiao (2001) proposed a model where both 
the dipole radiation and the longitudinal current outflow due to an 
unipolar generator take place. But the value of a they needed to fit 
the observed values of n did not agree with observationally deter- 
mined values of a. So even a better model is needed to understand 
pulsar magnetosheric properties better which might give even fur- 
ther different values of Bs than those obtained in this work. But 
the beauty of this work is that it reveals even orders of magnitude 
difference in the values of Bs from its canonical values within the 
simple dipole model. 
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5 CONCLUSION 

We have seen that depending upon EoS, Al and a, the value of Bs 
can differ significantly from the value quoted in the literature. The 
value of Bs also determines the value of the magnetic field through- 
out the pulsar magnetosphere by the relation B{r) oc Bs/r* where 
r is the radial distance from the star surface. So correct understand- 
ing of the dense matter physics will help us to understand pulsar 
magnetosphere related phenomena better. 

On the other hand, future technology might enable us to de- 
termine neutron star magnetic field strength {Bs) independently 
from the spectral analysis (like synchrotron absorption lines. Zee- 
man splitting etc). The values of a, M, P and P can be determined 
form timing analysis. Simultaneous determination of all these five 
quantities will provide an alternative way to constrain the dense 
matter EoS. 
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